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Abstract: In this paper we will present a mathematical description and analysis of 
a fractional-order regulated system in the state space and the state-space controller 
design based on placing the closed-loop poles on the complex plane. Presented are 
the results of simulations and stability investigation of this system. 
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1. Introduction 

Research of the fractional-order derivatives continued in the last decades not 
only as a mathematical branch [ft], % |I[ , etc., but in many applied domains, e.g. 
||, H and in control theory too Q - 16|, etc. In works j?], ||, ^| the first generalizations 
of analysis methods for fractional-order control systems were made (s-plane, frequency 
response, etc.). The following of the above-mentioned works were oriented to the 
methods of fractional-order system parameters identification, methods of fractional- 
order controllers synthesis, methods of stability analysis, methods of control of chaotic 
fractional-order systems, and so on. 
In works @ 

was presented a state space model described in vector and matrix 
relations expressing the fractional-order derivatives 

x (Q) (i) = A x(t) +B u(t), 

y(t) = C x(t), t > 0. (1) 

This description is convenient only for simple models of systems with only one 
fractional-order derivation, or for spatial type of systems. In work JT^ we proposed a 
state space model of the linear time- invariant one dimensional system which expresses 
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the first derivatives in the state space equations and has the classical state space 
interpretation for the fractional-order system too. On the right side of these equations 
we can then transfer more than one fractional-order derivatives of the state space vari- 
ables. A disadvantage of this expression is that in time domain we cannot express the 
state space equations in vector and matrix relations as in previous description. But 
we can do this in s-plane. 

This contribution deals with a mathematical description and analysis of a 
fractional-order regulated system in the state space and the state-space controller 
design based on placing the closed-loop poles on the complex plane. 

2. Definition of the system 

For the definition of the control system we consider a simple unity feedback control 
system illustrated in Fig.l, where G s (s) denotes the transfer function of the controlled 
system and G r (s) is the controller transfer function, both integer- or fractional-order. 
The differential equation of the above closed regulation system for the transfer func- 
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Figure 1: Unity- feedback control system 

tion of the controlled system G s (s) = l/(a 2 s a + a\s^ + ao) and the fractional-order 
PD S controller G r (s) = K + T d s s has the form 

a 2 y (a) (t) + a x yW (t) + T d (t) + (a +K)y(t) = K w(t) + T d (t) (2) 
where a, (3, 8 are generally real numbers and ao, a±, a^, K, T d are arbitrary constants. 

3. Fractional-order control system with PD S controller in the state space 

Consider the system described by differential equation (||) . After its modification 
and substitution of state space variables x(t) = Xi(t),x(t) — x 2 (t),x(t) = X2{t) we 
can derive the following state space model equivalent to model (||) 
x\(t) = x 2 (t) , 

Mt) = _^±J£. x f-<*) {t) _ Ti x 2 1+s - a \t) - ^ 4 1+/? - Q) (t) + - w ^- a \t) , (3) 

a 2 a 2 a 2 a 2 

y(t) = K xi(t)+T d 4 tf_1) (t) , t>0. 

Of course, we can make other alternative state space models for the same system. 
For the approximation of the fractional-order derivatives on the right-hand side of 
equations @ we can take the relation (ji|) from e.g. ||, [h]] 

N(t) 

y^(t)^h- a J2 b Mt~jh), (4) 

3=0 
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where L is "memory length", h is time step of calculation, N(t) = min j , j , 

[z] is the integer part of z, bj = (— 1) J ("), and ("J is binomial coefficient. To calculate 
bj it is convenient to use the following recurrent relation 



1 + ol . , 
bo = l, bj = (l — ) 



(5) 



After the above-mentioned fractional-order derivatives discretisation (^) and discreti- 
sation of the first derivatives on the left-hand side of equations (0) we obtained the 
simple Euler methods for solving the state space model 

x%,k+i = Xi,k + hx 2 ,k , 

k k 



X2,k+1 = X 2 ,k 



h - 



a 2 



t>jXl,k-j h > CjX 2 ,k-j 

do — 

3=0 Z j=0 



3=0 
k 

y k = K x\,k + T d h {1 ~ 5) ^2 ejx 2 ,k- 

3=0 



A-^X^fc-i + 1 ^ E 6 3^-3 



(6) 



3=0 



fc > 0. 



From these equations we can compute state trajectories of the fractional-order control 
system described above. 



4. Design of the PD S controller in the state space 

After Laplace transformation we can write state space equations (^) in vector and 
matrix relations in s-plane 



pX(a) = A(a) X(s) + B(a) W(a) 

y(s) = c(s) x(s), t > o 

where matrix A(s) and vectors B(s), C(s) are 
1 



(7) 



A( S ) = 



ao+-ftT „2-q 
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From equations (Q) we can derive the overall transfer function 

Y(s) _ C adj(sl — A) B 



W(s) 



det(s\ - A) 



(8) 



The characteristic equation of the closed-loop system is determined by solving the 
determinant in the denominator of the transfer function, which is a fractional-order 
polynomial in s 

a 2 s a +a 1 s l3 +T d s s + {a +K) = (9) 
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Solving this equation we can find the poles of the closed-loop system with known 
parameters, or design the PD S controller parameters for desired control system poles. 

Consider the system described by differential equation (||) with system coefficients 
a2 = 0.8, ai = 0.5, ao = I, a = 2.2,(3 = 0.9. The task is to determine the controller 
parameters K, Td and i5 for desired control system poles s\,2 = — l±6i and the steady 
state error less than 4%. Parameter K = 24 can be computed from the equation based 
on the steady state error e ss (oo) of the closed-loop control system K = (100/e ss — l)ao- 
If we introduce the desired system poles to the characteristic equation ([)]) we 
obtain system of two nonlinear equations from which we can easy derive two equations 
8 = arctan(2, 9839)^^ and T d 2 37 5 = 645, 2174 for the calculation of the controller 
parameters Td = 6, 9407 and 6 = 0, 71859. In Fig. 2 and Fig. 3 are depicted the unit 
- step responses of the classical numerical solution |l(J and the numerical solution of 
the state space model (Jfy. The obtained state trajectory represents stable focal point 
for the above-mentioned coefficients of the system. 




Figure 2: Unit step response Figure 3: State trajectory 

Assuming the integer-order PD controller we obtain a system of two linear equa- 
tions for direct calculation of the controller parameters K = 36, 0854 and T d = 4, 0141. 

If we increased in case of fractional-order PD S controller the requirement on steady 
state error less than 2% we would obtain the following controller parameters K = 49, 
T d = -79, 74427 and 6 = -0, 55194. That means, we obtained the controller with 
weak integrator. But the characteristic equation of such a system has one additional 
pole S3 = 1,98 and it follows from the stability analysis that such closed-loop system 
is unstable. We can verify this fact also with frequency methods for stability investi- 
gation of the fractional-order system jl^, [l3j . In this case we can change the desired 
control system poles Si 2, or we have to proceed in a different way, e.g. as follows. 

5. Fractional-order control system with PI X controller in the state space 

The differential equation of the closed regulation system for the transfer function 
of the controlled system G s (s) = l/(a,2S a + ais 13 + ao) and the fractional-order PI X 
controller G r (s) = K + TiS~ x has the form 

a 2 y {a+ ^ (t) + a, (t) + (a + K) (t) + T t y(t) = K (t) + T t w{t) (10) 

where a, [3, A are generally real numbers and ao, a\, K, Tj are arbitrary constants. 
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The state space model equivalent to model ( |l0| ) has the folloving form 
x\(t) — -x 2 {t) + w, 
X2(t) = +x 3 (t) , 

33(f) = +H xt a - X) (t) - xt a \t) - ^ X [ 1+l3 - a \t) + * W { t a \t) , (11) 

0,2 0-2 a 2 a 2 

y{t) = x 2 (t), t>0. 

where x \ (t) = z(t) is the actuating error signal and ^(i) = Vi(t), x$(t) — v-i{t) are the 
state variables of the controlled system (Fig. 4). As above, we can derive the following 




Figure 4: Block diagram representing the state and output equations 

characteristic equation of the closed-loop system 

a 2 s a+x + a 1 s l3+x + (a + K)s x + Ti = (12) 

Assume the same controlled system as in the previous section. If we introduce the 
desired system poles si, S2, S3 to the characteristic equation ( |l2|) we obtain a system 
of three nonlinear equations from which we can numcricaly calculate the controller 
parameters K, Ti and A. 

6. Conclusion 

We have presented a mathematical description of a fractional-order control system 
in the state space and the state-space controller design based on placing the closed- 
loop poles on the complex plane. In the design of a fractional-order PD S regulator 
it is necessary to pay attention to the question of stability of the control system. In 
the case of negative values of 5 the order of the system increases. In addition to the 
required poles, the system can develop new poles that can render the system unstable. 
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